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Abstract — We consider an improved rotational mechanism of the explosion of a collapsing supernova. 
We show that this mechanism leads to two-stage collapse with a phase difference of ^5 h. Based on this 
model, we attempt a new interpretation of the events in underground neutrino detectors on February 23, 
1987, related to the supernova SN 1987A. 
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INTRODUCTION 

In effect, the idea of two-stage gravitational col- 
lapse has long been a subject of discussion, par- 
ticularly in the case of the resumption (of the sec- 
ond stage) of the collapse in a neutron star with its 
transformation into a black hole (see, e.g., Imshen- 
nik and Nadyozhin 1988). However, if the rotation 
effects are taken into account, then such a two-stage 
collapse of an iron— oxygen— carbon (Fe— O— C) core 
acquires a different, more specific content that we 
have called the rotational mechanism of the explosion 
of a collapsing supernova (SN). The reason is that, 
in our opinion, the above rotation effects make it 
possible to solve the crucial problem of the trans- 
formation of the collapse into an explosion for high- 
mass and collapsing supernovae (all types of SN 
except the type-la thermonuclear SN). An extensive 
series of studies has been carried out since 1992 in 
connection with the famous SN 1987A (Imshennik 
and Nadyozhin 1992; Imshennik 1992; Aksenov and 
Imshennik 1994; Imshennik and Popov 1994; Ak- 
senov et al. 1995; Imshennik 1995, 1996; Imshen- 
nik and Blinnikov 1996). Below, we consider the so- 
called improved rotational mechanism of explosion. 
The possibility of this mechanism was first men- 
tioned by Imshennik and Popov (1994) in connec- 
tion with the reception of two neutrino signals from 
SN 1987A (see Section 1) separated by a relatively 
long time interval of 4 h 44 min (=4.7 h) (Dadykin 
et al. 1987; Aglietta et al. 1987) at t UT = 2 h 52 min 
(February 23, 1987) and t UT = 7 h 36 min (Hirata 
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et al. 1987; Bionta et al. 1987). These observations 
of the neutrino signals were carefully analyzed by 
Dadykin et al. (1989). The above improvement of the 
explosion scenario for the rotating Fe— O— C core of 
a high-mass star stems not only from the necessity 
of explaining the observations of the neutrino signals 
from SN 1987A, but also from the intrinsic logic of 
the development of this scenario — it arose on the road 
to overcoming theoretical difficulties. 

The neutrino spectra were obtained from theoret- 
ical estimates. These spectra are based on the hy- 
drodynamic calculations of a quasi-one-dimensional 
model for the formation of a rotating collapsar (Im- 
shennik and Nadyozhin 1977, 1992) and on the hy- 
pothesis of bulk neutrino radiation from a rotating 
collapsar (Imshennik and Nadyozhin 1972; Ivanova 
et al. 1969a) with the almost total dominance of 
electron neutrinos in the neutrino radiation (see Sec- 
tion 2). The derived spectra include the dimensionless 
chemical potential of electrons, cp, which is considered 
here as the only free parameter of these spectra. 

In Section 3, we analyze the observations of 
the neutrino signal on the LSD detector at t]jj = 
2 h 52 min (February 23, 1987) by using the previ- 
ously obtained neutrino spectra that additionally take 
into account the effects of self-absorption inside a 
rotating collapsar. We develop the hypothesis about 
the interaction of electron neutrinos with the nuclei 
of iron that is actually present in the LSD detector 
in large quantities. The products of these interaction 
reactions, mainly in the form of gamma-ray photons 
and electrons, are detected in a liquid scintillator 
with a photomultiplier. The detection efficiency is 
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estimated by the Monte-Carlo method. We show that 
the observational data are consistent with the theory 
described in Sections 1 and 2 over a wide <p range. 



1. THEORETICAL ANALYSIS OF THE TIME 
INTERVAL BETWEEN THE TWO STAGES 
OF COLLAPSE FOR SN 1987A 

Thus, on the threshold of gravitational collapse, 
the Fe— O— C stellar core has a given (from calcula- 
tions of the evolution of high-mass stars with a total 
mass on the main sequence M ms > 1OM0) mass Mt 
and a total angular momentum Jo, which are, obvi- 
ously, conserved during the collapse of this core into a 
rotating collapsar. Aksenov et al. (1995) numerically 
constructed a large family of such two-dimensional, 
axisymmetric collapsars as a function of the parame- 
ters Mt and Jo. Dong and Shapiro (1995) proved the 
high probability of collapsars falling into the region 
of dynamical instability that is specified by the stan- 
dard criterion f5 = £ ro t/|£gr a v| > 0.27(Tassoul 1978). 
The quantities £ mi and |£ gra v| denote the total rota- 
tional and total gravitational energies, respectively. 
Note that during collapse with the conservation of 
total angular momentum Jo and local specific angular 
momentum, the energy £ ro t greatly increases com- 
pared to |£g rav |, which, of course, is also an increasing 
quantity. This instability grows with a characteristic 
hydrodynamic time and typically leads to the breakup 
of the collapsar into pieces, in the simplest case, into a 
binary of neutron stars ( N S ); almost all of the angular 
momentum can transform orbital angular momen- 
tum, J or b < Jo. However, A J = Jo — J or b becomes 
the spin angular momentum of the NS themselves 
essentially in the more massive component of this bi- 
nary. In other words, the NS binary is formed through 
the hydrodynamic fragmentation of a rotating collap- 
sar. Imshennik (1992) showed that, under the addi- 
tional assumption of a circular orbit of the binary and 
for given Mt and J or b, all of the orbital parameters, in- 
cluding the orbital radius ao and velocity vq (in terms 
of the reduced binary mass MiM^jM^, according to 
the Kepler law) can be determined: 
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where the free parameter < 5o < 1/2 appears. Be- 
low, Mi is assumed to be the mass of the less mas- 
sive NS. Remarkably, the evolution of the binary is 
determined by only one factor — gravitational radia- 
tion. The latter is unique from an astrophysical point 
of view, because it pertains to the evolution of such 
binaries, but in the presence of a low-mass NS whose 



mass is much lower than the mass of the more mas- 
sive NS at the very outset; this mass is also small 
in absolute terms — compared to M© (see below). 
In the point-mass approximation, the gravitational 
radiation (Peters and Mathews 1963; Landau and 
Lifshitz 1973) and the conservative mass transfer are 
described by a simple differential equation (Imshennik 
and Popov 2002): 
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where, according to the common property of NS (de- 
generate stars), there is mass transfer with a decrease 
in the mass of the low-mass NS (d5/dt < 0). The 
derivative d5/dt can be explicitly calculated in the 
Roche approximation (the Roche lobe and potential), 
as was shown in detail in the paper mentioned above, 
which continues the classic works in this field pio- 
neered by Paczynski and Sienkiewicz (1972). Equa- 
tion (2) then describes the evolution under the action 
of two factors; the second factor is mass transfer 
that causes the radius a to increase, in contrast to 
the first factor that describes the approach of the 
components. However, the second factor comes into 
play only after the low-mass NS fills its Roche lobe. 
Before this time, the NS only approach each other at 
a constant parameter 5 = 5q due to the gravitational 
radiation that carries away not only the energy of the 
NS binary, but also its orbital angular momentum. 
For d5/dt = 0, i.e., at 5 = 5o, Eq. (2) has a simple 
analytical solution from which i grav , the time of the 
closest approach of the components with a constant 
arbitrary parameter 5o, can be obtained formally up to 
the radius a = (Imshennik and Popov 1994): 
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where jo = J or b/8.81 x 10 49 erg s and m t = 
Mt/2M0; this choice of scales for the dimensionless 
quantities jo and m\ corresponds to the typical 
conditions of a Fe— O— C stellar core on the threshold 
of its collapse (as applied to SN 1987A). It should 
be noted that the arbitrariness in specifying the initial 
rotation in the stellar core is severely restricted by the 
hypothesis that this rotation is rigid; this is related 
to the action of convection inherent in the structure 
of high-mass stars, particularly at the final stages of 
their thermonuclear evolution. For this reason, the 
rotational energy in the initial conditions is actually 
negligible compared to the gravitational energy, i.e., 
j3<l, and the stellar structure is virtually spherically 
symmetric. 

Table 1 presents the values of 5q that follow from 
relation (3), as long as m t = 0.9 (M t = 1.8M©) and, 
most importantly, as applied to SN 1987A, t grav = 
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Table 1 



./orb, ergs 


8.81 x 10 49 


6.17 x 10 49 


3.17 x 10 49 


jo 


1.00 


0.700 


0.360 


Jac, ergs 


1.72 x 10 49 


1.45 x 10 49 


1.25 x 10 49 




0.195 


0.235 


0.394 


Mi,M0 


0.37 


0.25 


0.13 


So 


0.206 


0.139 


0.0722 


At gr , S 


-0.04 


-0.2 


-10 


At ac , s 


-0.9 


-1.0 


-1.2 



Note: M t = 1.8M0, J = 8.81 x 10 49 erg s, t grav = 4.7 h. 



4.7 h = 16 920 s are given. In these calculations, we 
vary jo, i.e., the fraction of the total angular momen- 
tum Jo transformed into the orbital angular momen- 
tum J or b. According to the data in Table 1, jo varies 
between 1 and 1/3. The next row gives the orbital 
angular momentum J ac at the time the low-mass NS 
fills its Roche lobe and the evolution of the binary with 
mass transfer begins. This evolutionary stage is called 
accretion for short; clearly, the solution of Eq. (2) 
with the second term d5/dt / on its right-hand side 
corresponds to this stage. We see that J ac decreased 
significantly at the previous stage of NS approach 
under the effect of gravitational radiation alone: by a 
factor from ~5 to ~3, depending on j . Finally, the 
next rows give the sought-for values of the mass Mi 
and 5q (its ratio to the total mass Mt), which also 
decrease by a factor of almost 3. It is easy to see that 
the following strong inequality holds over the entire jo 
range: Mi <c M2 (at least, M2 is larger than Mi by a 
factor of 4 for the first column). 

Thus, our identification of the time interval be- 
tween the two neutrino signals for SN 1987A with the 
time of approach of the components of the putative 
NS binary due to gravitational radiation alone (based 
on relation (3)) fits into the theory of the rotational 
mechanism of supernovae explosions for the follow- 
ing reason: The minimum NS mass Mi = 0.13MQ 
still exceeds the lower mass limit for stationary NS, 
Mi min = 0.095M© (Blinnikov et al. 1990; Imshen- 
nik 1992; Aksenov et al. 1995), when the star ex- 
plosively disintegrates with its transformation into an 
iron gas. The derived strong inequality 5q <C 1 quali- 
tatively agrees with the first three-dimensional hydro- 
dynamic calculations of the fragmentation of a rotat- 



ing collapsar (Houser et al. 1994; Aksenov 1999), in 
which a ~0.1M© mass ejection actually emerges. 1 

Nevertheless, an analysis of the physical effects 
disregarded in relation (3) is required to justify the 
results presented in Table 1 more reliably. To this 
end, the last two rows of Table 1 give Ai gr and 
Ai ac that have the following meaning: The values of 
these quantities should be compared with the £ grav = 
16 920 s specified above. Our estimates separately 
are smaller than the latter value at least by a factor 
of 2000 and are of no importance within the accuracy 
of the sought parameter 5q. The physical meaning 
of Atg r is that, strictly speaking, from the time tg rav in 
relation (3), we should subtract the time interval when 
the radius of the NS binary a = oro, where or , the 
critical radius for filling the Roche lobe for a low- mass 
NS, is(Paczynski 1971) 

a R0 = 2.16i? N si<5 - - 33 , (4) 

where R^si is the radius of a NS with mass Mi. 2 
Thus, it follows from (4) that a R0 > -Rnsi at <5o <C 1. 
It is easy to show that At gr = t grav ( Jac/Jorb) 8 - The 
values of Ai gr are given in Table 1 ; they are very small. 
Imshennik and Popov (1998) analyzed in detail the 
mixed type of evolution described by the complete 
equation (2) with mass transfer (!) until S decreased 
to its critical value of (5 cr = Mi min /M t = 0.053, i.e., 
until the explosion time of the low-mass NS. By 
definition, mass transfer is possible if the condition for 
the immersion of a low-mass NS in its Roche lobe 
is satisfied, i.e., a < a#, which includes the current 
critical radius for filling the Roche lobe, 

a R = 2.16i? NS <r a33 , (4') 

obtained by a natural generalization of relation (4), 
depending on the parameter 5 < 5q and the ra- 
dius i?Ns of an NS with mass Mi = 5M{. This 
inequality is actually satisfied throughout the mixed 
evolutionary stage of the NS binary, which can also 

1 Interestingly, if a rotating collapsar fragmented into pieces of 
equal mass (So = 1/2) at the same values of m t = 1.8 and 
jo = 1.0, the time of approach of the components (t grav ~ 
400 s), according to formula (3), would be many times 
shorter than the time specified by the observations of the 
neutrino signals (tgrav = 16 920 s)! 

2 As the NS radius Rns, we may take interpolation formulas 
derived for numerical calculations of the radii of cold NS with 
the inclusion of low-mass NS, for examples, from Yaranows- 
ki and Kxolak (1992). However, it should be remembered 
that, in our case, still very young and, hence, relatively hot 
NS whose radii are generally larger than the radii of cold NS 
are members of the binary. This implies that their Roche 
lobes are filled earlier and that their radii <zro are larger, but 
the related changed corrections At gr and At ac are still negli- 
gible compared to the sought-for time t gra v from relation (3). 
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be proved from physical considerations (Imshennik 
and Popov 1996). 

Ai ac follow from these calculations, because the 
calculations yield a time dependence At(6) within the 
range 5q > S > <5 cr with At ac = At(S cr ). In contrast 
to At gr , this quantity should be added to i grav , because 
the two effects partly compensate each other. There- 
fore, their total effect, which is at a maximum in the 
last row of Table 1 , is 8.8 s, which, we repeat, is negli- 
gible for the sought-for parameter 5q. In short, quan- 
titatively including the second evolutionary stage of 
the components of the NS binary turned out to be 
completely unnecessary. 

Nevertheless, it remains to determine how impor- 
tant the assumptions about the total mass of the 
Fe—O—C stellar core, My (which was taken above to 
be 1.8M0), about the total angular momentum Jo 
(which was taken above to be 8.81 x 10 49 erg s), and, 
finally, about a zero orbital eccentricity eo = are. At 
a fixed tgrav, obvious dependences follow from (3) for 
<5 < 1: (1 - S )S ~ <5 oc j§- 89 m t -1 - 7 . In Table 1, the 
dependence So oc Jq 89 , of course, holds. Since the de- 
crease in initial Jo as well as its increase are severely 
restricted by the condition of the dynamical rotational 
instability itself (Aksenov et al. 1995), the change in 
Mi oc J^ 89 oc J 9 - 89 is small. The same is also true 
for the dependence Mi oc M t -1 ' 7 that follows from the 
previous dependence of the dimensionless quantities, 
but for a different reason: the masses of Fe—O—C nu- 
clei are limited to the standard range 1.2M0 < Mt < 
2M0, as predicted by the stellar evolution theory. 
Note that, in this case, M t = 1.8M0 in Table 1 is 
close to the upper limit of the range, so the possi- 
bility of its decrease in stellar cores before collapse 
definitely does not bring the corresponding values 
of So outside the inequality <5o > S cr - We can reach the 
important conclusion that all the possible changes in 
parameters M t and J do not qualitatively change the 
results that were formulated above based on the data 
of Table 1. 

Next, let us consider not a circular but an eccentric 
orbit that may well result from the hydrodynamic 
fragmentation of a collapsar, i.e., with an initial ec- 
centricity e / (0 < eo < 1). The problem of the 
evolution of an eccentric orbit for an NS binary was 
completely solved by Imshennik and Popov (1994). 
Although its solution (also analytical at an arbitrary 
value of eo) is more complex than that in the limiting 
case of a circular orbit with e = 0, the quantity i grav 
of interest can be expressed elegantly (see the pio- 
neering paper by Peters 1964) at the same values 
of M t and J as those in (3): 



tgrav 
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where the new (compared to (3)) dimensionless func- 
tion r(eo) is given in the form of an easily calculated 
integral with a certain factor (both are functions of the 
parameter eo): 

, v 48 1 
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This function r(eo) > 1 differs only slightly from unity 
(t(0) = 1) as long as e < 0.5 and increases steeply 
at eo — > 1, r(eo) oc 1.81(1 — e^) -1 / 2 . Our analysis 
shows that the limiting (singular) case of eo = 1 
itself cannot be considered, because the necessary 
condition for the orbits being quasi-stationary is not 
satisfied. For this reason, this case has been excluded 
in the above eo range. The qualitative result obtained 
by Imshennik and Popov (1994) may be considered 
to be the establishment of such a rapid decrease 
in eccentricity e (de/dt < for any e and So) due 
to gravitational radiation that e = e\ < 0.1 by the 
time the low-mass NS fills its Roche lobe, i.e., the 
orbit is almost indistinguishable from a circular orbit. 
To be more precise, this inequality also depends on 
the parameters my and j and on the radius R^si 
(Imshennik and Popov 1994): 



*< 0044(f) 



19/12 
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13.5 km 



19/12 



(7) 



As our analysis indicates, the influence of these pa- 
rameters and the radius is quantitatively small, so 
the previous eccentricity estimate is definitely justified 
by formula (7). This circumstance seems quite fortu- 
nate, because the entire theory of mass transfer was 
constructed precisely for circular orbits. It would be 
instructive to estimate the influence of a finite initial 
eccentricity eo on Sq. We take jo = 1 and my = 0.9 
in formula (5) and substitute eo = 0.9, so the func- 
tion r(eo) = 3; hence we obtain So = 0.241 at i grav = 
4.7 h, which should be compared with Sq = 0.206 
(the first column in Table 1 ). The following general 
conclusion can be reached: the influence of a finite 
initial eccentricity, eo / 0, is negligible and does not 
lead (like the influence of my) to any violation of the 
condition So > <5 cr . The latter inequality ensures the 
existence of a low-mass ejection (So <C 1) in the form 
ofanNS. 

One of the important results of the analytical 
model for the evolution of a close NS binary is the 
conclusion that, in the solution of Eq. (2) with mass 
transfer, the approach of the components very soon 
gives way to their recession from one another, and the 
dependence on the initial arbitrary parameter Sq of the 
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binary virtually disappears. We are talking here about 
the sought-for dependence of 5 on time At during 
the mixed evolutionary stage under the action of 
both factors: gravitational radiation and mass transfer 
(accretion) (see Fig. 6 in the paper by Imshennik 
and Popov (1998) cited above). This circumstance 
is formally similar to what has been said above 
about the influence of the other initial parameter eo 
on the solution. Thus, during the evolution under 
consideration, the influence of both emerged arbitrary 
parameters, 5q and eo, essentially vanishes, and the 
evolution of a close NS binary finishes irrespective 
of the choice of these parameters. However, a third 
arbitrary parameter, A J = Jo — J or t>, appeared in the 
course of our theoretical analysis. This parameter 
may prove to be crucial: Will the evolution end with 
the merger or explosion of a low-mass NS that 
has overfilled its Roche lobe? As Imshennik and 
Popov (1998) showed, the concept of a (circular) 
orbit and the point approximation for the gravitational 
interaction of both NSs become meaningless at the 
very end of the evolution. 3 A numerical solution 
of the three-dimensional hydrodynamic problem is 
required. This solution is extremely complex if the 
effects of gravitational radiation, neutrino radiation, 
and the nonideal equation of state for these NS are 
included in the problem. Of course, this solution is of 
relevant interest, but we have to restrict our selves 
to reasonable (?) physical estimates before it can 
be implemented after the overcoming of fantastic 
difficulties. Until now, it has been assumed that 
the more massive component with mass Mi = (1 — 
<5o)M t does not collaps for the second time during 
the fragmentation of the rotating collapsar due to 
the residual total angular momentum AJo. The first 
neutrino signal was produced only by the primary 
collapse of the Fe— O— C stellar core that led to the 
formation of this collapsar before its fragmentation. 
However, even during the evolution, once the orbital 
radius a = a^o was reached, the high-mass NS 
gradually (in a time At ac ~ 1 s) increased its mass up 
to Mi = (1 — S cr )Mi. However, the centrifugal forces 
attributable to the angular momentum A J could still 
reliably prevent its secondary collapse. Here, it is 
important to note that this angular momentum AJ 
can even exceed J ac : for example, for the second 
column of Table 1, J ac = 1-45 x 10 49 erg s, while 
the corresponding AJ = 2.64 x 10 49 erg s. Using 
relations ( 1 ), but for the final state of a binary with 

Quantitatively, this last period of the evolution of a close NS 
binary is very short in terms of 8; more specifically, it arises 
at 8 = 8* = 0.082, which exceeds the above critical value 
of (5 cr = 0.053 only slightly. The parameter 8* is virtually 
independent of the initial value of 8o, but depends on other 
parameters, M t and J or b (Imshennik and Popov 1998). 



the orbital angular momentum J ac given above, we 
can easily find that the parameters of this orbit are 
a f = 2.22 x 10 7 cm and v t = 3.28 x 10 9 cm s _1 , so 
the angular frequency of the binary is uj\ = Vf/af = 
1.48 x 10 2 s -1 . On the other hand, the angular 
frequency of a high-mass NS with the above value 
of M 2f = 1.70M© is fi f ~ 10 5 s" 1 for rigid rotation 
(overestimate) or Of ~ 10 3 s~ x for a typical differen- 
tial rotation law (underestimate). 

In short, the strong inequality Of » u\ holds. 
However, the corotation condition is valid in close 
binaries. According to this condition, the orbital 
angular momentum must be much larger than the 
spin angular momenta, because the frequencies of 
these rotations, u\ and Of, are equalized by tidal 
forces. 

Strictly speaking, this condition takes effect even 
at the very beginning of the mixed evolutionary stage, 
i.e., since the time the low-mass NS fills its Roche 
lobe. Thus, a corollary of the corotation condition is 
the transformation of spin angular momentum AJ 
into the orbital angular momentum of the binary. In 
this way, we also remove this main obstacle to the 
secondary collapse of a high-mass NS. The following 
legitimate question arises: Why did not the corota- 
tion condition prevent the appearance of a significant 
residual angular momentum A J during the fragmen- 
tation of the rotating collapsar? A completely justified 
answer to this question will probably be given after 
the construction of a three-dimensional (!) hydro- 
dynamic model with a consistent allowance for the 
same physical factors that we mentioned when dis- 
cussing the end of the evolution of a NS binary, but we 
can make a qualitative estimate. To be more precise, 
we can draw attention to the characteristic times of 
the hydrodynamic fragmentation, which was called 
hydrodynamic not by chance. These characteristic 
times are only 10~ 3 s, while the characteristic times 
of the mixed evolution for a NS binary, as has been 
repeatedly said above, are much larger, ~1 s. There 
is yet another factor in favor of the appearance of the 
residual angular momentum under consideration in a 
high-mass NS, as will be seen below. 

Thus, a high-mass NS collapses in complete ac- 
cordance with the standard model without rotation 
(Nadyozhin 1977a, 1977b, 1978). This collapse is 
accompanied by a standard neutrino signal, because 
during its evolution, the binary has virtually gotten 
rid of the bulk of its total angular momentum due to 
gravitational radiation. However, the remnant of the 
spin angular momentum was redistributed, accord- 
ing to the above assessment, to the orbital angular 
momentum, essentially transferring to the low-mass 
NS (see also below). It is important to recall that the 
integrated parameters of the standard neutrino signal 
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are in satisfactory agreement with these parameters of 
the second neutrino signal from SN 1987A (Imshen- 
nik and Nadyozhin 1988; Blinnikov et al. 1988). 

2. ESTIMATING THE (ELECTRON) 
NEUTRINO SPECTRA FOR SN 1987A 

So far, we have been able to consistently ana- 
lyze the collapse of a stellar Fe— O— C-core with ini- 
tial rotation only in the one-dimensional, spherically 
symmetric statement of the problem by averaging the 
centrifugal force over the meridional angle (the only 
change of the problem then concerns the equation of 
motion from the complete system of four equations 
written in Lagrangian coordinates). The assumption 
that the local specific angular momentum j is con- 
served played a key role in this great simplification 
of the problem. In this case, the centrifugal force F r 
was unambiguously expressed in terms of the initial 
rigid-rotation parameters ujq = const and ro(m) and 
the Eulerian radius r = r(m,t), with F r oc r~ 3 . Ob- 
viously, the total angular momentum was then also 
automatically conserved in this case, 

Jo = ^ J jdm = J rl(m)dm, 
o o 

with the factor 2/3 being attributable precisely to the 
angular averaging of the centrifugal force: 

2wgr$(m) 



3 r 



3 r 3 (m, t) 



(8) 



The first numerical calculations of this kind were 
performed by Imshennik and Nadyozhin (1977). The 
equations of state, the processes of material neu- 
tronization, and the description of neutrino processes 
did not differ in any way from their representation 
in the standard model. These hydrodynamic models, 
which were called quasi-one-dimensional because of 
the described allowance for the rotation effects, were 
computed until the formation of a hydrostatic equi- 
librium configuration that was called above a rotating 
collapsar. 

Another study of the quasi-one-dimensional model 
with detailed analysis and testing of the results was 
carried out later (Imshennik and Nadyozhin 1992). 
In this paper, we focus on the discussion of the 
neutrino radiation parameters that were limited in the 
papers by Imshennik and Nadyozhin (1977, 1992) 
only to data on a light curve with a total energy 
6 V y = 3.3 x 10 52 erg (by the end of the computation 
t = 2.9 s) and to parameters of the neutrinosphere. 
Apart from the marked difference between the neu- 
trino light curves for the calculations of the stan- 
dard model (cjq = 0) and the quasi-one-dimensional 



model with uj = 0.86ft («o = (GMt/R^) 1 / 2 is the 
velocity at which the centrifugal and gravitational 
forces are equal on the surface of the stellar core 
near the equator), 4 an enormous difference appears 
between the neutrino optical depths indicated by 
several numbers in these light curves (see Fig. 6 from 
Imshennik and Nadyozhin 1992). Whereas r > 100 
for the standard model starting from the maximum of 
the light curve, r is typically several units (except the 
middle part of the light curve where it rises to 32) 
in the model with rotation. This implies that the 
approximation of radiative heat conduction used in 
the calculations is near its validity boundary for the 
quasi-one-dimensional model. In addition, the region 
inside the neutrinosphere encloses less than half of M t 
by mass (see Fig. 1 from the cited paper). By the end 
of the computation, this mass decreases by several 
more times (~5) (see Fig. 5 from the cited paper). 
We must also take into account the ambiguity in 
determining the neutrino optical depth itself (formally, 
four different determinations of it are possible!), in 
contrast to the photon optical depth. Nevertheless, 
an optically opaque region near the center of the 
collapsing stellar core since the time t = 0.5 s was 
introduced in the cited paper (see Fig. 1). The main 
objection to this introduction is as follows: The rotat- 
ing collapsar in the initial two-dimensional axisym- 
metric configuration is a highly flattened structure 
(Aksenov et al. 1995) the polar radius of which is 
several times smaller than its equatorial radius. This 
flattening naturally entails a decrease in r at least by 
the same number of times. In addition, it is clear that 
the rotating collapsar will be a hydrostatic equilibrium 
configuration only in an axisymmetric geometry. If 
the criterion for dynamical instability given above 
and obtained quantitatively, f3 = 0.42 > 0.27 = /3 cr 
(Imshennik and Nadyozhin 1992), is satisfied, in- 
stability relative to the third (azimuthal) coordinate 
will rapidly (on the characteristic hydrodynamic time 
scale) transform it into a dumbbell-like configuration. 
The central region of the collapsar with the lowest 
specific angular momenta will be located in the bar 
of this dumbbell, which is very thin compared to 
the end balls. Here, it should be emphasized that, 
strictly speaking, the rotating collapsar cannot be 
formed as an equilibrium structure if the criterion 
[5 > f3 cr is satisfied. It has the right to exist only for 
our idealization of axial symmetry and, of course, in 



4 Since the adopted initial angular velocity uio is of great im- 
portance, we give its numerical value specified in the calcula- 
tions of the quasi-one-dimensional model: Qo = 1.78 s -1 (at 
M t = 2M0 (!), Ro = 4.38 x 10 8 cm), and the correspond- 
ing value of loo = 0.86fio = 1.53 s -1 (the substitution of the 
coefficient 0.86 for 0.80 in the succeeding paper (Imshennik 
and Nadyozhin 1992) is among the corrections to the pre- 
ceding paper (Imshennik and Nadyozhin 1977) 
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terms of the quasi-one-dimensional model. Thus, 
the former center of the collapsing stellar core will 
certainly be seen not only toward the poles, but also 
on the sides of the bar, into which it falls during 
the nonlinear growth of dumbbell barlike (m = 2) 
instability (Tassoul 1978; Aksenov 1996). 

The above qualitative discussion gives us the right 
to advance the following hypothesis: the central re- 
gion of the quasi-one-dimensional model with the 
formation of a rotating collapsar may be considered 
to be completely transparent for intrinsic neutrino 
radiation. Concurrently, it seems quite justifiable to 
retain the thermodynamic parameters in this region 
obtained in the calculations by Imshennik and Na- 
dyozhin (1977, 1992) as moderately sensitive to this 
change of the status of neutrino radiation (from sur- 
face to bulk), because the optical depths are so small 
in these calculations. In any case, retaining such 
quantities as the central density and temperature (p c 
and T c ) we can be sure of their consistency, because 
they satisfy three of the four equations of the quasi- 
one-dimensional model (together with a violation of 
the entropy equation). Concurrently, we retain the 
total energy of the neutrino radiation, £*- = 2.7£ u t = 
8.9 x 10 52 erg. The numerical coefficient 2.7 is in- 
troduced by analogy with the standard model (Nady- 
ozhin 1977a, 1977b, 1978), where the corresponding 
value of £ V y = 1.9 x 10 53 erg at time t = 2.4 s (seethe 
table from Imshennik and Nadyozhin 1992) trans- 
forms into £* v = 5.3 x 10 53 erg through justified esti- 
mates of the prolonged neutrino cooling stage of a hot 
NS (Imshennik and Nadyozhin 1988). The presented 
approach may be fraught with the overestimation of 
the bulk radiation parameters, particularly the total 
time of this process for an unstable rotating collap- 
sar. 5 The interpretation of the observed time interval 
between the two neutrino signals from SN 1987A 
justified in the preceding section on the basis of the 
previously advanced rotational mechanism of the ex- 
plosions of collapsing supernovae may serve as an 
indirect justification of the suggested approach. 

Thus, we take the total energy carried away by 
neutrino radiation with allowance made for the cool- 
ing of the rotating collapsar (by the time ~ 6 s), 

£*. = 8.9 x 10 52 erg. (9) 

At this time, the material is so cold that the neu- 
trons degenerate (actually, this effect should have 
been taken into account at the end of the numerical 
calculation of the quasi-one-dimensional model, t = 
2.9 s). 



5 We are grateful to D.K. Nadyozhin for a critical discussion 
in which he, in particular, drew special attention to the above 
circumstance. 



Let us next consider the spectral properties of 
the neutrino radiation by assuming the bulk radi- 
ation mechanism due to the non-one-dimensional 
hydrodynamic collapse noted above. We restrict our 
analysis to the main reaction of the modified URCA- 
process e~ + p — > n + v e , in which electron neutri- 
nos v e are generated. If the material contains only free 
nucleons (see the estimates below), as was specified 
in the quasi-one-dimensional model, then, based on 
the Kirchhoff law and using the expression for the 
v e mean free path (Ivanova et al. 1969a; Imshennik 
and Nadyozhin 1972), we obtain the spectral specific 
radiation power: 

q u _ 4ttB u 1 In 2 1 
p p m {ft) np 1 + 6 



\m e c 2 J l + exp(^-v?) 

where 6 = N n /N p is the degree of neutronization of 
the material, ip is the dimensionless chemical poten- 
tial of the electron gas (see below), e v is the neutrino 
energy in erg, and the spectral power is given per unit 
interval of this energy. Below, we will be interested in 
the dependence of (q u /p) on energy e„, in which an 
allowance for the Pauli principle in the initial v e mean 
free path plays a crucial role: it is responsible for the 
exponential factor in (10). Note also that relation (10) 
gives the spectral power of the neutrino radiation B v 
per unit solid angle. By definition, the neutrino radi- 
ation in the reaction under consideration is isotropic 
in each (Lagrangian) particle of the material, so B v 
in (10) is also isotropic. The energy dependence of 
(q„/ p) follows from (10): 

-ocf^y (in 

p \m e c z J 

x 77 r ~- —? r = <l>{ x i ¥>) > 

1 + exp (f- - </>) l + exp(x-p) 

where x = e u /kT, p = n e /kT, and \x e is the chemical 
potential of the electrons ([p e ] = erg)- The function 
4> = 4>{x, p) with a specified parameter p has only one 
maximum with x max > p and the following asymp- 
totics: ~x 5 (for x — > 0) and ~ x 5 e~ x (for x — > oo). 
The quantity x max = x max (p) can be easily calculated 
and is given in Table 2. Curiously, at ip = 10, the exact 
solution for x max = p = 10, while in the remaining 
cases of specified p>, x max is given with four significant 
figures. We probably could have also calculated the 
v e energy averaged over spectrum (11), but it will 
suffice to understand that it is close to x max , and 
exceeds it slightly. These important quantities in the 
v e spectrum appear to have the following main prop- 
erty: x m ax > 5 for p < 5; x max ~ p> for p > 5). The 
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inequality x max > 5 always holds for the lower limit 
of the maximum of the spectrum. 

It would be appropriate to apply a correction for 
the self-absorption of neutrinos of sufficiently high 
energies in the layers surrounding the center of the 
rotating collapsar to the v e spectrum from (11). Let 
the mean optical depth of these layers be (t v ) = k, 
where k is an arbitrary number that may be close to 
unity (k < 1) even if the non-one-dimensional effects 
are taken into account. Of course, it is smaller than 
the optical depths obtained by Imshennik and Nady- 
ozhin (1992) (see above). The value of (t u ) can be 
determined by the following estimate together with its 
spectral value: 



Table 2 



<r„) = R±- 



f B u (l/Qde u 



f B u de u 
o 



R 



(12) 



where R is an effective radius (to simplify the estima- 
tion of ( 12), we disregard the weak dependence of the 
thermodynamic parameters in the surrounding layers 
on the current radius r). Thus, the mean (t„) in (12) 
is weighted in radiation intensity B v from (10). Let us 
write the expression for the mean free path which 
will be needed below (recall that it is uniquely related 
to ( 1 0) for B v by the Kirchhoff law): 



ra 1 + 6 1 

cr 9 p 



m P c 



1 + exp(e„/(fcr) - tp) 

exp(e„/(fcT) - <p) ' 
(13) 



whence we find an explicit expression for the mean (t u ): 



(r„) = R- 



kT 



too 1 + \ m e c 
The dimensionless function \I/(</?) is given by 



*(¥>)• (14) 



1 + exp(x — <p) 



dx 



(15) 



x 



1 + exp(x — <p) 



dx 



-i 



This expression can be easily tabulated by numerical 
integration. Table 2 gives the values of ^(ip). In con- 
clusion, we can easily determine the dimensionless 
energy y = e v /(kT) that with l' v from (13) separates 
the spectrum into two parts: r u < kiix < y and t„ > 
k if x > y. To this end, we equate the spectral depth t v 
to the mean (t„) from (14). Substituting Eq. (13) 
for l' v and Eq. (14) for (t u ) into the equality t„ = (t u ) 








2.5 


5.0 


7.5 


10.0 


20 




5.033 


5.303 


6.327 


8.008 


10.00 


18.97 


*(v>) 


42.30 


44.29 


50.14 


59.70 


71.39 


125.9 


y(<p) 


6.509 


6.750 


7.396 


8.757 


10.57 


19.32 



yields a transcendental equation for the sought-for 
quantity y: 

y 2 = *(<p)[e X p(<p-y) + i\. (16) 

The numerical solution of Eq. ( 16) is also presented in 
Table 2. We see that y exceeds x max only slightly for all 
values of <p, implying that the v e spectrum at k ~ 1 is 
cut off immediately after its maximum. For k < 1, the 
cutoff boundary is shifted slightly to higher energies. 
Note that the true value of k can be determined in 
future hydrodynamic calculations, but its closeness 
to unity follows from the paper by Imshennik and 
Nadyozhin(1992). 

Below, we apply these simple estimates to the 
thermodynamic parameters of the quasi-one-dimen- 
sional model for a rotating collapsar. Unfortunately, 
only some of the characteristic quantities, and only 
in the form of p c and T c , the central density and 
temperature, but not the parameter ip c , the central 
chemical potential of the electrons from (10), are at 
our disposal. 6 To calculate the corresponding value 
of ip c , we will do the following. We will use the elec- 
trical neutrality condition for the material, which is 
naturally satisfied in the hydrodynamic calculations of 
collapse. It is greatly simplified in the ultrarelativistic 
case for a Fermi— Dirac electron gas: 



8tt fkT\ . o 



mo(l + 0)' 



(17) 



We can determine <p of interest from this condition 
(in the form of a cubic equation) by using the given 
parameters p and T. We know these parameters at 
the stellar center, p c and T c , for the time t = 2.9 s 
at which the computation ends (Imshennik and Na- 
dyozhin 1992; see the table there) and for the ear- 
lier time t = 0.0 s at which neutrino opacity sets in 
(Imshennik and Nadyozhin 1977). Two pairs of these 
parameters are presented in Table 3. However, as 
we see from condition (17), the parameter 6 should 



6 In both their publications of the calculations of the quasi- 
one-dimensional model, Imshennik and Nadyozhin (1977, 
1992) focused on the rotation effects and the possibility of 
the transformation of collapse into an explosion and provided 
only minimum information concerning the integrated quan- 
tities of the light curve for neutrino radiation and its total 
energy. 
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Table 3 



Pc 


= 2.6 x 10 14 gcrrr 3 , 
T c = 6.2 x 10 10 K 
(kT c = 5.34 MeV) 


Pc 


= 1.70 x 10 13 gem- 3 , 
T c = 5.89 x 10 10 K 
[kT c = 5.074 MeV) 


e c 


12.6 


100 


1000 


e c 


12.6 


(fie 


25.7 


13.0 


5.63 




10.6 



also be specified to solve the cubic equation. In Ta- 
ble 3, it is varied for the first pair of p c and T c over 
a wide range; the lowest value from the previous 
case will suffice for the second pair. In reality, the 
hydrodynamic calculations yield values of the material 
neutronization parameter 6 that do not exceed 100. 
Therefore, we conclude from the data of Table 3 that 
ip c ~ 10 may serve as a rough estimate. The corre- 
sponding energies, the mean energy and the energy at 
the maximum of the spectrum, in ( 1 1 ) are then found 
to be equal (see Table 2): 

(e v ) ~ £„ m ax ~ kT c ip c (18) 
~ (53.4-50.7) MeV ~ 50 MeV. 

Estimate (18) should be improved further, but the 
unusual hardness of the neutrino spectrum compared 
to the standard model of collapse without rotation 
is beyond question. We emphasize that the steep 
fall attributable to the v e self-absorption considered 
above (see Table 2) takes place immediately after the 
maximum of the spectrum. 

Next, let us first estimate the total number of v e 
required for the almost complete neutronization of the 
material of a rotating collapsar (6 < 100) indepen- 
dent of the mean energy of the v e spectrum: 

N v = 1.8M — ( ^] = 1.0 x 10 57 . 
m \A Fe J 

On the other hand, according to the previous estimate 
of N v , the total energy £ v of these neutrinos with 
(e v ) = 50 MeV is 

£ v = N v {e v ) = 8.0 x 10 52 erg, 

a value that is only slightly lower than £ v p from (9). 
To be more precise, about 10% of the energy remains 
for antineutrinos, v e : 

£ = £ v0 -£„ = 0.9 x 10 52 erg. 

Note that the necessary condition £ v „ > £ v is still 
satisfied, with nothing remaining for muon and taon 
neutrinos {£ VjJl vn = £ V tvt = 0). Meanwhile, the large 
deficit of v e was clear from our estimate of the elec- 
tron chemical potential. According to Imshennik and 
Nadyozhin (1972), the latter is related to the neutrino 



chemical potential by a simple relation for the ther- 
modynamic equilibrium conditions of the v e {y e ) and 
e~(e + ) degenerate gases: 

ip = if — In 0. 

Thus, for ip = 10, we obtain ip = 4.6—7.7, depending 
on = 100—10. This implies that the number of i> e is 
strongly suppressed, being proportional to exp(— ip), 
compared to the number of v e . 

It may seem that we use the parameters T c and p c 
noncritically, with the density of the material (on the 
order of nuclear density!) playing a major role in the 
achievement of high characteristic values of ip c ~ 
10. However, the density p c was carefully analyzed 
by Imshennik and Nadyozhin (1992), including its 
dependence on the initial values of the rotation pa- 
rameter loq, using the derived analytical formula that 
relates this density to the relative centrifugal force 
near the center of the rotating collapsar. This analy- 
sis inspired confidence that initial rigid rotation with 
a sufficiently large total angular momentum, Jo ~ 
10 50 erg s, does not prevent the emergence of a high 
density, p c ~ 10 14 g cm~ 3 , at the center of the rotat- 
ing collapsar. 

One would think that a significant admixture of 
iron nuclei etc. may be retained in the material under 
these conditions of relatively cold collapse (for a com- 
parison of the central temperatures of the calcula- 
tion under consideration and the standard model, see 
Imshennik and Nadyozhin (1992)). This admixture, of 
course, would immediately made condition (17), from 
which the chemical potential ip was estimated, more 
complex. The iron mass fraction X Fe can be estimated 
using the paper by Imshennik and Nadyozhin ( 1 965), 
where a universally accepted equation of state for 
the material was derived. Thus, using formula (29) 
from the cited paper for the specified parameters from 
Table 3, we find that X Fe /X n ~ 1.3 x 10" 2 ; i.e., for 
X n ~ 1, we actually obtain a negligible iron mass 
fraction, X Fe ~ 10~ 2 . Note that a correction on the 
order of unity to condition (17) then arises, but it also 
implies an increase in the sought-for parameter <p. 

In conclusion, we may mention one subtle effect 
called the redistribution of specific angular momen- 
tum. In the quasi-one-dimensional model, this redis- 
tribution, of course, is constant, but it can change 
in principle during the fragmentation of a rotating 
collapsar. This would even be necessary at So ~ 0.5; 
i.e., when the collapsar is divided into equal pieces. 
What physical factors could provide this redistribu- 
tion? First, we may point out the impact of three- 
dimensional tidal forces on any shear viscosity (for 
estimates of the neutrino viscosity, see Imshennik 
and Nadyozhin (1992)) and, finally, on the magnetic 
field. Of course, these factors should be included in 
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Table 4 



Model 


fi.erg 


£2, erg 


£3, erg 


ep e ,MeV 


£„ e ,MeV 


^, T ,MeV 


T, s 


I 


(3-14) x 10 53 


(0.5-2.3) x 10 53 


1Q 52 


12.6 


10.5 




20 




II 








10 


8 


25 


5 



Note. £\ is the total burst energy transformed into neutrinos of all types; £2 is the total energy carried away by Vi, where Ui — 
v e , v e , Vp, v^, v T , v T \ £3 is the total energy carried away by v e during the neutronization of the star in a time ~3 x 10~ 2 s; es e , 
e Ve , e Vll T are the spectrum-averaged energies of v e , v e , and v^ T , respectively; T is the duration of the neutrino burst. 



future hydrodynamic models. At present, however, we 
may dispense with them if we are dealing with frag- 
mentation into two pieces with a large difference in 
their masses (So <C 1). The initial rigid-rotation law 
then provides a high orbital angular momentum of a 
low-mass ejection almost without any redistribution 
of specific angular momentum. In our opinion, this 
circumstance is another argument for fragmentation 
with low values of the parameter <5n (see Table 1) 
suggested in Section 2. 

It would probably make sense to discuss the 
main conclusion of this section concerning the hard 
spectrum of electron neutrinos in the absence of other 
types of neutrinos from the viewpoint of possible 
future studies. We repeat that this conclusion requires 
a new series of quasi-hydrodynamic calculations of 
collapse in which, first, the bulk nature of the neutrino 
radiation in the generalized URCA process with 
corrections for self-absorption (Ivanova et al. 1969a) 
would be postulated, and, second, the more general 
kinetic equation for neutronization to which the 
complete system of kinetic equations reduces in 
nuclear statistical equilibrium (Imshennik and Na- 
dyozhin 1982) would be used in the calculations of 
the degree of neutronization of the material instead of 
the rough approximation of kinetic equilibrium for the 
/3-processes. We may take the risk of asserting that 
such calculations will generally lead to a decrease 
in the central temperature of the rotating collap- 
sar, but with a simultaneous increase in its central 
density. According to relation (18), such changes in 
thermodynamic parameters will entail approximate 
conservation of or, most likely, a moderate decrease in 
the characteristic energies of the v e spectrum, which 
will probably be in better agreement with the LSD 
observations of the neutrino signal from SN 1987A 
than the very hard spectrum with (e v ) ~ e^max — 
50 MeV in ( 18) considered above. 



3. A MODEL FOR A ROTATING COLLAPSAR 

AND A POSSIBLE INTERPRETATION OF 
THE EXPERIMENTAL RESULTS OBTAINED 
WITH NEUTRINO DETECTORS AT THE 
EXPLOSION TIME OF SN 1987A ON 
FEBRUARY 23, 1987 

Let us consider how the various detectors operated 
during the explosion of SN 1987A could record the 
neutrino signals in terms of the model for a rotating 
collapsar that reduces to the following: 

( 1 ) Two neutrino bursts separated by a time t grav ~ 
5 h must exist. 

( 2 ) The neutrino flux during the first burst consists 
of electron neutrinos (u e ) with a total energy £ Ue = 
£„ e = 8.9 x 10 52 erg; the neutrino energy spectrum 
4>(£ Ue ,(p) ( 1 1 ) is hard and asymmetric with mean en- 
ergies in the range 25—50 MeV (Fig. 1 ); the duration 
of the neutrino radiation is ig n ~2.9— 6 s. 

(3) The second neutrino burst corresponds to the 
theory of standard collapse. 

To compare the detector responses to these two 
neutrino bursts, we first give basic parameters of the 
neutrino fluxes obtained in the standard model and 
consider the neutrino detection methods. 

The experimental searches for neutrino signals 
from collapsing stars began with the paper by Zel- 
dovich and Guseinov (1965), who showed that gravi- 
tational collapse is accompanied by an intense short 
pulse of neutrino radiation. The search for collapse 
by detecting such a signal was first suggested by 
Domogatsky and Zatsepin (1965). The role of neu- 
trinos in stellar collapses was considered by Arnett 
(1966), Ivanova et al. (1969b), Imshennik and Na- 
dyozhin (1982), Nadyozhin and Otroshchenko (1982) 
[model I], Bowers and Wilson (1982), Wilson et al. 
(1986) [model II] and Bruenn (1987). 

Parameters of the neutrino fluxes during the col- 
lapse of nonmagnetic, nonrotating, spherically sym- 
metric stars were obtained by the above authors and 
are given in Table 4. 

It follows from the theory of standard collapse that 
the total energy carried away by neutrinos of all types 
(u e , i> e , Vp, P M , u T , v T ) corresponds to ~0.1 of the 
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Fig. 1. Neutrino energy spectrum §{e v , ip) (in arbitrary units) versus e v in MeV (a) and 4>*(x, tp) (in arbitrary units) versus 
x — e v /kT (b); tp — ^e/kT is the dimensionless chemical potential of the electron gas. 



stellar core mass and is equally divided between these 
six components. 

Until now, Cherenkov (H2O) and scintillation 
(C n H2„) detectors, which are capable of recording 
mainly i> e , have been used in searching for and de- 
tecting neutrino radiation. This choice is natural and 
is related to the large cross section for the interaction 
of £> e with protons: 

v e +p-> e + + n, (19) 

a 0eP 9.3£f+ x 1(T 44 cm 2 , E e+ » 0.5 MeV. 

Moreover, reaction (19) has a distinctive signature, 
and, as was first shown by Chudakov et al. (1973), 



a proton may be used as a neutron catcher with the 
subsequent formation of deuterium (d) and the emis- 
sion of a 7-ray photon with a time r ss 180—200 11m: 

n + + £ 7 = 2.2MeV. (20) 

Electron antineutrinos v e can be detected by search- 
ing for a pair of pulses separated by the capture time. 
The first and second pulses are attributable to the 
detection of e + and 7-ray photons, respectively. 

The z/ e -scattering reactions have much smaller 
cross sections, but they allow the neutrino arrival 
direction to be determined in Cherenkov detectors: 
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v e + e~ = u e + e~, a Ue e = 9Ae Ve x 10" 45 cm 2 , e Ve > 0.5 MeV, (21a) 

z> e + e~ = P e + e~, cr/> ee = 3.9e„ e x 10~ 45 cm 2 , e„ e > 0.5 MeV, (21b) 

\ Vi + e - = Ui + e - = i.6e x 10 -45 cm 2 , > 0.5 MeV, . , 01 , 

S a* o i = h,t. (21c) 
+ e" = P e + e", cr^ e = 1.3e^ x 10~ 45 cm 2 , e^>0.5MeV, 

In scintillation detectors, neutrinos with energies above the threshold (-Ethr) a ' s ° interact with carbon: 

Vi + 12 C -> 12 C* + z/j, ^thr = 15.1 MeV, i = e,n,r, (22a) 
L^ 12 C + 7 , £ 7 = 15.1 MeV, 

zv e + 12 C -> 12 N + e", E thr = 17.34 MeV, (22b) 
L+. 12 C + e + + i/ e , t= 15.9 ms, 

P e + 12 C -> 12 B + e+, E t hr = 14.4 MeV, (22c) 
L^ 12 C + e~ + u e , t = 29.3 ms, 



where E ihr is the threshold energy of the reaction, and 
r is the mean lifetime of the 12 B and 12 N isotopes. 

The cross sections for reactions (22a)— (22c) were 
calculated by Donnely (1973), Fukugita et al. (1988), 
Mintz et al. (1989), Kolb et al. (1994), and Engel 
et al. (1996) and measured in the Los Alamos E225 
(Athanassopoulos et al. 1990), LSND (Allen 1997), 
and KARMEN experiments (see, e.g., Mashuv 1998). 

Reaction (22a), which yields a monochromatic line 
at 15.1 MeV, allows one to measure the fluxes of 
and v T with energies twice the v e energy (Ryazhskaya 
and Ryasny 1992) for standard collapse and the to- 
tal neutrino flux with energy E v > 15.1 MeV for any 
model of collapse. Reactions (22b) and (22c) in the 
standard scenario yield a smaller number of events 
without oscillations. With oscillations, the number 
of events increases and becomes comparable to the 
effect produced by reaction (22a) at equal neutri- 
nosphere temperatures for neutrinos and antineutri- 
nos, T Pe (Aglietta et al. 2002). 

The cross section for the interaction of neutrinos 
with 16 nuclei in Cherenkov detectors a(u 16 0) for 
e v < 25 MeV is low compared to the cross sec- 
tion a{v 12 C) (Bugaev et al. 1979). 

In Cherenkov and scintillation detectors, the 
emission from e + , e~ , and photons is recorded by 
a photomultiplier. A neutrino burst is identified by 
the appearance of a series of pulses in the range of 
amplitudes from the detector energy threshold (E t h r ) 
to 50 MeV in a time from several seconds to several 
tens of seconds, depending on the model of collapse. 

Four detectors were operating during the explo- 
sion of SN 1987A: two scintillation (in USSR and 
Italy) and two Cherenkov (in USA and Japan). Pa- 
rameters of the detectors are given in Table 5. 



The effects from reactions (19) and (21a), (21c) 
expected in the model of standard stellar collapse in 
the Large Magellanic Cloud are shown in Table 6. We 
see that reaction (19) gives the largest contribution. 
It should be noted that all of the detectors except 
the LSD could record only e + ; i.e., there was no 
signature of the event by which reaction (19) could 
be unambiguously identified. 

Features of the Construction of Scintillation 
Detectors 

Since we explore the possibility of detecting elec- 
tron neutrinos by the nuclei of the detector materials, 
we will briefly consider the features of the construc- 
tion of scintillation detectors, because some of them 
have not been reported previously. 

The Baksan Underground Scintillation Telescope 
(BUST) (Alexeyev et al. 1979) is a cube with sides 
14 m square. There are two horizontal planes inside 
the cube. 400 scintillation counters 0.7 x 0.7 x 0.3 m 3 
in size are located on each of the cube's planes 
(six outer and two inner). Each counter is watched 
by one photomultiplier. The spacing between the 
horizontal planes separated by an absorber is about 
3.6 m; the thickness of the absorber is ~170 g cm~ 2 , 
~20 g cm~ 2 of which is made up of iron. The scintil- 
lator is made of white spirit (Voevodskii et al. 1970); 
the molecular composition is C n H 2 „, n = 10. The 
detection threshold of energy release in the module 
is ^10 MeV. The working volume of the detector 
consists of counters located on three horizontal 
planes: the two inner and one lower; the scintillator 
mass is 132 1. In analyzing the results of February 23, 
1987, the authors also used data from some of the 
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Depth of 
water 
equivalent, m 




Detection 
threshold, 
MeV 


Detection efficiency 


Background 
pulse 
frequency 
m.s- 1 ** 


Detector 


Working mass, t 
Material 


e + spectrum 
of reaction 
v e p — ► e + n 
(19) 


e~ spectrum 
of reaction 

v l e" — > z^e" 
(21a, 21c)* 


BUST, USSR 


850 


130(200) C„H 2n 
160 Fe 


10 


0.6 


0.15(0.54) 


0.013(0.033) 


LSD, 

USSR-Italy 


5200 


90 C„H 2 „ 

200 Fe 


5-7 


0.9 


0.4(0.7) 


0.01 


KH, 

Japan— USA 


2700 


2140 H 2 


7-14 


0.7 


0.17(0.54) 


0.022 


1MB, USA 


1570 


5000 H 2 


20-50 


0.1 


0.02(0.18) 


3.5 x 10~ 6 



* The detection efficiencies of the electron spectrum produced in the reactions v^ipn^) + eT — > ^^(^.t) + e~ (21c) are given in 
parentheses. 

** The background is given in the energy range E thr — 50 MeV; for the Cherenkov detectors, the background is given for the recording 
of internal events. 



outer detectors, which increased working mass to 
200 t. 

The LSD detector (Badino et al. 1984; Dadykin 
1979) operated under Mont Blanc at a depth of 
5200 m.w.e. consists of nine modules located on three 
floors with an area of 6.4 x 7.4 m 2 . The LSD height 
was 4.5 m. The module is an iron container with 
an area of 6.4 x 2.14 m 2 , a height of 1.5 m, and a 
wall thickness of 2 cm, with eight cells separated 
by 2-cm-thick iron sheets in which 1 x 1.5 x 1 m 3 
scintillation counters are located. In fact, LSD is an 
iron scintillation detector. To reduce the influence of 
radioactivity from the surrounding rocks, the detector 
is shielded by steel plates. The total iron mass is about 
200 t. The mass of the scintillator is made up of white 



Table 6 



Detector 




K e - (21a) + (21c) 


K e - (21c) 


LSD 


1.5 


0.043 


0.024 


BUST 


2 


0.052 


0.036 


KH 


17 


0.53 


0.36 


1MB 


6 


0.4 


0.35 



spirit (Voevodskii al. 1970) is 90 1. Each counter is 
watched by three photomultipliers (FEU-49B). The 
high sensitivity and low background of the detector 
allow both e + and n particles in reaction (19) to be 
detected. 

The parameters of the Kamiokande II and 1MB 
detectors are well known (see Hirata et al. 1987; 
Bionta etal. 1987). 

Experimental Data 

To understand how the recorded events during the 
collapse of SN 1987A on February 23, 1987, fit into 
the scenario of a rotating collapsar, we present them 
in the chronological sequence of observation: 

(1) February 23, 1987, 2:52:36.79 UT: LSD 
recorded a cluster of five pulses in real time; the 
estimated probability of this being mimicked by 
background radiation is low (less than once in three 
years). No such event had been observed in two years 
of LSD operation (1985-1987) (Aglietta et al. 1987; 
Dadykin et al. 1987, 1988). At present, we may add 
that the detector recorded no similar signal over its 
entire operation time until 1999. However, it is worth 
noting that the background conditions in 1988 were 
improved by an additional shield. 

(2) February 23, 1987: information about the event 
measured by the LSD was transmitted to the group 
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February 23, 1987 



1 

1 1 


3 5 

i i i 


7 

i i 


9 11 

i i i i 


Optical observations 






| m v = 12 m 




m v = 6 m | 


Geograf 


2:52:35.4 


LSD 5 

■ 


2:52:36.8 
43.8 


2 


7:36:00 
19 


™ (4 ) = 


2:52:34 
44 


' 2 1 


7:35:35 
47 


1MB 




8 a 


7:35:41 
47 


BUST , 


2:52:34 


6 B 


7:36:06 
21 



Fig. 2. Time sequence of the events recorded by various 
detectors on February 23, 1987. 

working with the GEOGRAV gravitational antenna 
(Rome). 

(3) February 23, 1987, 10 : 40 UT: the first optical 
observation is made of SN 1987A in the Large Mag- 
ellanic Cloud, ~50 kpc from Earth, and is reported 
only the next day (IAU Circ. 1987). 

It is worth emphasizing that the LSD event was 
observed and reported to Amaldi's group before in- 
formation about the supernovae explosion was ob- 
tained. The data from the other detectors were ana- 
lyzed later. Figure 2 shows the time sequence of the 
events recorded by various detectors on February 23, 
1987. We see that there are two groups of events 
concentrated near the LSD time (2:52:36) and the 
1MB and KN time (7:35:35). All the events have been 
extensively discussed for the last several years. 

Interestingly, strong correlations between the 
pulses from gravitational antennas and underground 
detectors were observed from 2:00 UT until 8:00 UT 
(Amaldi et al. 1987; Pizzella 1989; Aglietta et al. 
1989, 1991), but the analysis of this fact is not the 
subject of this paper. Here, we wish to analyze the 
event recorded under Mont Blanc in terms of the 
model of a rotating collapsar. The events measured by 
underground detectors near 7:36 UT are interpreted 
by most physicists as the detection of antineutrino 
radiation, while attempts to explain the LSD data at 
2:52 UT in a similar way run into difficulties. 

Table 7 presents the parameters of the LSD pulses. 
The third column gives energies of the measured 
pulses corresponding to the two muon calibrations 
made before and after February 23, 1987 (each in a 
time interval of three months) and separated by two 
months. We see that the energy is determined with an 
accuracy of 20—25%. 




e v , MeV 



Fig. 3. Comparison of the reduced total cross sections 
with the neutrino cross section on a free neutron for the 
reaction v e + (A, Z) — > e~ + (A, Z + 1)* (Bugaev et al. 
1979). 

Difficulties in Interpreting the Effect Measured 
by LSD at 2:52 UT in the Case 
of Antineutrino Detection 

As was shown by Dadykin et al. (1989), a to- 
tal neutrino radiation energy of E v = 6£„ e ps 1.2 x 
10 55 erg is required to explain the LSD effect in terms 
of the detection of antineutrino radiation if the con- 
flicts with the results of other underground detectors 
have been eliminated. This energy is more than an 
order of magnitude higher than the binding energy of 
a neutron star with a baryon mass of about 2M Q . 

In addition, only one of the five measured pulses in 
the cluster was accompanied by a neutron-like pulse 
that was offset from the trigger signal in the detector 
by 278 ixm and that had an energy of 1.4 MeV. As- 
suming the detection of five antineutrinos, one would 
expect the detection of five positrons accompanied, on 



Table 7 



Event no. 


Time, UT ± 2 ms 


Energy, MeV 


1 


2:52:36.79 


7-6.2 


2 


2:52:40.65 


8-5.8 


3 


2:52:41.01 


11-7.8 


4 


2:52:42.70 


7-7.0 


5 


2:52:43.80 


9-6.8 


1 


7:36:00.54 


8 


2 


7:36:18.88 


9 
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Table 8 





a = 


1.27 x 10" 


An 9 
-40 cm 2 


K.e 


= 31.84 MeV, 










E~y = 


1.82 MeV, ny. J2 n E 1 = 1.72 MeV 




a = 


6.41 x 10" 


- 41 cm 2 


Ek, s - 


= 30.84 MeV, E 1 = 1 MeV, 










Ey = 


1.82 MeV, n 7 : J2 n E i = 1 - 72 MeV 


GT 


a = 


1.05 x 10" 


-40 cm 2 


E K,e- 


= 27.84 MeV, E 1 = 4 MeV, 










E~y = 


1.82 MaB, nr. £„ ^7 = L72 MeV 


GT 


a = 


1.27 x 10" 


-40 cm 2 


Ek,c- 


= 24.84 MeV, E 1 = 7 MeV, 










E-y = 


1.82 MeV, nr. £„ £ 7 = 1.72 MeV 



* The Fermi level. 
** The Gamow-Teller resonance. 
*** The electron kinetic energy. 



average, by two neutrons. In this case, the probability 
of measuring one pulse attributable to the neutron 
capture by hydrogen and offset from the positron sig- 
nal by 1.5t, where r is the neutron lifetime, is less 
than 5%. 

The above arguments make the explanation of the 
LSD effect in terms of antineutrino detection implau- 
sible. 



A Possible Explanation of the LSD Effect in the 
Model of a Rotating Collapsar 

As was mentioned above, the collapsar emits elec- 
tron neutrinos (u e ) with a total energy of £v e * ~ 
8.9 x 10 52 erg, the spectrum shown in Fig. 1, and 
mean energies of 30—40 MeV for ~2.9— 6 s. These 
neutrinos can be recorded by the detector nuclei via 
the reactions 

L e + (A,Z)^e- + (A,Z + l) 

\i/ c + (A,Z)-e- + (A,Z + l)*, 

v e + (A,Z)^i/ e + (A,Z)*. (23b) 

The detectors operated on February 23, 1987, con- 
tained either oxygen, mainly 16 (KH and 1MB), or 
carbon, mainly 12 C, and iron 56 Fe (LSD, BUST). 
It follows from the paper by Bugaev et al. (1979) 
that the reduced cross section a Uen = a Ue A/N for iron 
at e v < 40 MeV exceeds ov e „ for oxygen by more 
than a factor of 20 (ov e „( 56 Fe) > 20a„ en ( 16 0)) (see 
Fig. 3). Thus, for these energies, the number of u e A 
interactions in the LSD (200 t of Fe) will be larger 
than that in the KH (1900 1 16 0). 

To answer the question of how the u e Fe inter- 
actions are detected in the LSD, we turn to Ta- 
ble 8. As an illustration, this table presents the partial 



cross sections calculated for the following reaction for 
e v = 40 MeV (Gaponov et al. 2003): 

i/e + iiFe->i?Co*+e- (24) 

The ground level of ^Fe— + , the ground level of 
27C0— 4 + . The difference between the binding ener- 
gies is E&Co] - E[§Fe] = 4.056 MeV. 

The threshold energy for reaction (24) is 8. 1 6 MeV. 
An electron can be produced with an energy from 
~3 1 .8 to 24.8 MeV, and its appearance is always 
accompanied by cascade 7-ray photons with a total 
energy of from 3.54 to 10.54 MeV. Recall that the 
critical energy (e) in iron (the electron energy at which 
the ionization losses are equal to the radiative losses) 
is 21 MeV. Thus, an electron with an energy E > e in 
iron on a thickness of d > 1 1 units (t is the radiation 
unit of length, 1 1 units = 13.9 g cm -2 (1.78 cm)) 
produces a small electromagnetic cascade. Calcu- 
lations indicate that, during the interaction of elec- 
tron neutrinos with e Ve = 40 MeV in a 2— 3-cm-thick 
iron layer located between two scintillation layers, 
many more 7-ray photons than electrons fall into the 
scintillator (Dedenko and Fedunin 2003). The mean 
energies of these particles are ~7— 9 MeV. The de- 
tection efficiency of z^ e Fe interactions (77) depends on 
the design of the detector and on the energy thresh- 
old: r\ ~ 75% for the inner part of the LSD (~90 t 
of Fe), r] ~ 35% for its outer part (~1 10 t of Fe), and 
r] < 15% for the BUST. We see that in the scenario 
under consideration, the mean energies recorded by a 
scintillation detector with an iron interlayer are close 
to the energies measured by the LSD on February 23, 
1987. The estimated effect of the detection of neutrino 
radiation in the first phase of the collapse of a rotating 
star by different detectors is presented in Table 9. 
We used the cross sections of reactions (23a) and 
(23b) for 56 Fe and 12 C and the cross sections of 
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Table 9 



Detector 


Detection 
threshold 


Estimated number of v e A interactions 


Estimated 
enect 


Experiment 






N 3 


N 4 


N 2 7] 


LSD 


5-7 


3.2 


5.7 


3.5 


4.9 


3.2 


5 


KH 


7-14 


0.9 


3.1 


1.2 


2.5 


2.7 


2* 


BUST 


10 


2.8 


5.2 






-1 





* See DeRujula( 1987a, 1987b). 
** Alexeyeve/a/. (1987). 



reaction (23a) for 16 (Bugaev et al. 1978; Gaponov 
et al. 2003; Fukugita et al. 1988; Haxton 1987). The 
estimates were obtained for monochromatic neutri- 
nos with e„ e = 30 MeV (JVi), e„ e = 40 MeV (JV 2 ) 
and for the electron neutrino spectrum (Fig. 1) at 
ip = 5 (N 3 ) and ^ = 7.5 (7V 4 ). 

We see that these estimates are consistent with the 
experimental data. If almost only electron neutrinos 
with a mean energy of ~35— 40 MeV were emitted in 
the first collapse, then the experimental data corre- 
spond to the scenario for the rotational mechanism of 
a supernova explosion. 

To be able to detect neutrino radiation from future 
stellar collapses (if the above scenario is realized), it is 
necessary to have detectors capable of recording not 
only v e , but also v e with high efficiency. Two detectors 
of this type exist: the LVD (1.1 kt of scintillator, 1 . 1 kt 
of Fe) and the SNO ( 1 kt of D 2 0). The Super K and 
Kamland detectors will be able to clearly see v e with 
e„ e > 40 MeV. 

CONCLUSIONS 

We have put forward and partly justified our inter- 
pretation of the signal recorded by the underground 
LSD detector at 2:52 UT on February 23, 1987, as 
the first detection of a neutrino burst from the col- 
lapse of SN 1987A. We proceeded from the previously 
suggested rotational mechanism of the explosions of 
collapsing supernovae and the idea of taking into 
account the interaction of electron neutrinos with the 
nuclei of iron whose presence in the LSD's construc- 
tion enormously increases the sensitivity to the first 
phase of collapse. A careful study of the correspond- 
ing nuclear reactions for the deneutronization and 
excitation of iron nuclides when electron neutrinos of 
sufficiently high energies, mainly from 20 to 50 MeV, 
interact with them led us to this conclusion. From the 
viewpoint of the rotational mechanism, it may be as- 
serted that a rotating iron stellar core collapses in two 
phases separated by a relatively longtime interval; the 
first phase of collapse was most likely detected by the 
LSD. According to this mechanism, the first phase 



is peculiar in that a rotating collapsar is formed in 
it with the emission of a very hard electron neutrino 
spectrum attributable to the reaction e~ + p — ► n + 
v e with an energy at the maximum of the spectrum 
up to 50 MeV in the almost complete absence of 
electron antineutrinos and other types of neutrinos 
(muon and taon). A more detailed study of this spec- 
trum based on numerical calculations (quasi-one- 
dimensional model) not only allows us to confirm the 
above properties but also to justify why other neutrino 
detectors have not recorded the first neutrino sig- 
nal. We additionally took into account other possible 
nuclear reactions of electron neutrinos with oxygen 
nuclides in the KH and 1MB Cherenkov detectors. 

In addition, according to the same mechanism, we 
can reliably interpret the time interval between the 
first and second neutrino signals. This time interval 
is mainly attributable to the well-known gravitational 
radiation generated during the fragmentation of the 
rotating collapsar into a binary of NS with greatly 
differing masses. The evolution of such a binary was 
theoretically analyzed in detail in Section 1. At this 
time, gravitational radiation carried away a signifi- 
cant fraction of the initial angular momentum of the 
iron core, which was undoubtedly conserved during 
the first phase of collapse. It may be asserted that 
allowance for the rotation effects of the collapsing iron 
core alone allows us to theoretically interpret the two 
successive neutrino signals from SN 1987A. Until 
now, the presence of a neutrino signal on the LSD 
detector has been ignored in the theory of this famous 
event, which is besides enigmatically correlated with 
the signals from two gravitational antennas: in Italy 
and the United States. We have not yet been able 
to reasonably interpret how the mentioned antennas 
responded to gravitational radiation during the frag- 
mentation of the rotating collapsar into a NS binary, 
which was theoretically estimated to be very modest. 

It remains to be added that the second neutrino 
signal recorded by the KH, 1MB, and BUST was as- 
sociated with the detection of electron antineutrinos 
predicted in the universally accepted standard model 
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for the secondary collapse of the high-mass NS in 
the putative binary once it had accreted the bulk of 
the mass from the low-mass NS and gotten rid of 
the remnants of its angular momentum. Theoretically, 
this signal corresponds to the standard hydrodynamic 
one-dimensional model of collapse without rotation 
with the formation of a neutrinosphere and with an 
equal energy distribution between all types of neutri- 
nos. Electron neutrinos could not be recorded on the 
LSD detector at this time, because these neutrinos 
must have much lower mean energies, about 15 MeV. 
The cross section for their interaction with iron nu- 
clides is too small. 

Of crucial importance is the fact that a low-mass 
NS simultaneously ceases to exist once it has reached 
a critical mass of about O.1M and has been de- 
stroyed by an explosion relatively far (several hundred 
km) from the high-mass NS. This transformation 
into an ejection of iron with the release of recombina- 
tion energy (~5 MeV/nucleon) and in the presence 
of high kinetic energy (~0.3 x 10 51 erg) leads to an 
explosion with quite a sufficient total energy of about 
10 51 erg, and with a directed symmetry (the direction 
of the initial orbital motion of the ejection); i.e., it 
solves the fundamental problem of the transformation 
of iron-core collapse into a supernova explosion with 
an observed total energy of about 10 51 erg. 

ACKNOWLEDGMENTS 

Undoubtedly, we would be unable to reliably inter- 
pret the five famous events recorded by the LSD with- 
out the great progress in calculating the cross sec- 
tions for the interaction of neutrinos with iron nu- 
clides, or to reliably justify the high detection effi- 
ciency of neutrino interactions in the LSD without 
calculating the development of low-energy electro- 
magnetic cascades. We wish to thank Yu.V. Gaponov, 
S.V. Semonov, L.G. Dedenko, and E.Yu. Fedunin, 
who made these calculations available to us. The 
strong stimulating influence of G.T Zatsepin, whom 
we wish to thank, on the development of the theory 
of collapse and supernova explosion and the experi- 
mental methods for detecting neutrino radiation from 
these natural phenomena, including SN 1987A, is 
hard to overestimate. We wish to thank L.B. Okun 
for his interest in the problem, D.K. Nadyozhin for 
his helpful critical discussions, D.V. Popov for his 
help in calculating the evolution of a NS binary, and 
V.A. Matveev for his invariable interest in this prob- 
lem. We are grateful to all members of the LSD 
collaboration, especially to C. Castagnoli, G. Cini- 
Castagnoli, V.L. Dadykin, V.G. Ryasny, O. Saave- 
dra, and W. Fulgione for their fruitful discussions; to 
V.B. Braginskii, G. Pizzella, and E. Coccia for their 
discussion of the possibility of recording gravitational 



radiation during SN 1987A, and to N.Yu. Agafono- 
va, N.A. Vulikh, and V.V. Kuznetsov for their help 
in preparing the paper for publication. This work 
was supported by the Russian Foundation for Basic 
Research (project nos. 03-02- 164 14a and 03-02- 
16436) and grants from scientific schools (projects 
NSh- 1787.2003.2 and NSh-1782.2003.2). 

REFERENCES 

1. M. Aglietta et al., Nuovo Cimento C 12,75(1989). 

2. M. Aglietta et at., Nuovo Cimento C 14, 171 (1991a). 

3. M. Aglietta et al., Nuovo Cimento B 106, 1257 
(1991b). 

4. M. Aglietta et al., Europhys. Lett. 3, 1321 (1987). 

5. M. Aglietta et al, Nucl. Phys. B (Proc. Suppl.) 1 10, 
410(2002). 

6. A. G. Aksenov, S. I. Blinnikov, and V. S. Imshennik, 
Astron. Zh. 72, 717 (1995) [Astron. Rep. 39, 638 
(1995)]. 

7. A. G. Aksenov and V. S. Imshennik, Pis'ma Astron. 
Zh. 20, 32 (1994) [Astron. Lett. 20, 24 (1994)]. 

8. A. G. Aksenov, Pis'ma Astron. Zh. 25, 163 (1999) 
[Astron. Lett. 25, 127(1999)]. 

9. A. G. Aksenov, Pis'ma Astron. Zh. 22, 706 (1996) 
[Astron. Lett. 22,634(1996)]. 

10. E. N. Alexeev et al., Pis'ma Zh. Eksp. Teor. Fiz. 45, 
461 (1987) [JETP Lett. 45, 589(1987)]. 

11. E.N. Alexeev et at., Proc. of the 16th ICRC.-Kyoto, 
Japan, 10,282(1979). 

12. R. C. Allen, Phys. Rev. Lett. 64, 1871 (1990). 

13. E. Amaldie/ 1 at., Europhys. Lett. 3, 1325(1987). 

14. W. D. Arnett, Can. J. Phys. 44, 2553(1966). 

15. C. Athanassopoulos et al., Phys. Rev. C 55, 2078 
(1997). 

16. J. Badino et at., Nuovo Cimento 7, 573 (1984). 

17. R. M. Bionta et at., Phys. Rev. Lett. 58, 1494(1987). 

18. S. I. Blinnikov, V. S. Imshennik, D. K. Nadyozhin, 
et al., Astron. Zh. 67, 1181 (1990) [Sov. Astron. 34, 
595(1990)]. 

19. S. I. Blinnikov, V. S. Imshennik, and D. K- Nadyozhin, 
Astrophys. Space Sci. 150,273(1988). 

20. R. Bowers and J. R. Wilson, Astrophys. J. 263, 366 
(1982). 

21. S. Bruenn, Phys. Rev. Lett. 59, 938(1987). 

22. E. V. Bugaev, G. S. Bisnovatyi-Kogan, et at., Nucl. 
Phys. A 324, 350 (1979). 

23. A. E. Chudakov, O. G. Ryajskaya, and G. T. Zatsepin, 
Proc. of 13th ICCR, Denver, 1973, Vol. 3, p. 2007. 

24. V. L. Dadykin, G. T. Zatsepin, and O. G. Ryazhskaya, 
Usp. Fiz. Nauk 158, 140 (1989) [Sov. Phys. Usp. 32, 
459(1989)]. 

25. V. L. Dadykin et al, Proc. of the 16th ICRC, Kyoto, 
Japan 10,285(1979). 

26. V. L. Dadykin et al., Pis'ma Astron. Zh. 14, 107 
(1988) [Sov. Astron. Lett. 14,44(1988)]. 

27. V. L. Dadykin et at., Pis'ma Zh. Eksp. Leor. Fiz. 45, 
464 (1987) [JELP Lett. 45, 593 (1987)]. 

28. L. G. Dedenko and E. Yu. Fedunin, Private Commu- 
nication (2003). 



A ROTATING COLLAPSAR 



18 



29. G. V. Domogatsky and G. T. Zatsepin, Proc. of 9th 
ICCR, London, 1965, Vol. 2, p. 1030. 

30. Dong Lai and S. L. Shapiro, Astrophys. J. 442, 259 

(1995) . 

31. T. W. Donnely, Phys. Lett. B 43B, 93 (1979). 

32. J. Engel et al., Phys. Rev. C 54, 2740 (1996). 

33. M. Fukugita, Y. Kohyama, and K. Kubodera, Phys. 
Lett. B 212, 139(1988). 

34. Yu. V. Gaponov, O. G. Ryazhskaya, and S. V. Se- 
menov, Yad. Fiz. (2003, in press). 

35. W. C. Haxton, Phys. Rev. D 36, 2283 (1987). 

36. K. Hirata et at., Phys. Rev. Lett. 58, 1490 (1987). 

37. J. L. Houser, J. M. Centerella, and S. C. Smith, Phys. 
Rev. Lett. 72, 1314(1994). 

38. IAUCirc. No. 4316(1987). 

39. V. S. Imshennik, Pis'ma Astron. Zh. 18, 489 (1992) 
[Sov. Astron. Lett. 18, 194(1992)]. 

40. V. S. Imshennik, Space Sci. Rev. 74, 325 (1995). 

41. V. S. Imshennik, Eruptive Stars, Ed. by A. G. Mase- 
vich (MID RF, Moscow, 1996), p. 7 [in Russian]. 

42. V. S. Imshennik and S.I. Blinnikov, IAU Coll. 145: 
Supernovae and Supernova Remnants, Ed. by 
R. McCray and Zhenru Wang (1996), p. 1 19. 

43. V. S. Imshennik and D. K. Nadyozhin, Astron. Zh. 42, 
1154(1965). 

44. V. S. Imshennik and D. K. Nadyozhin, Zh. Eksp. Teor. 
Fiz. 63, 1548(1972). 

45. V. S. Imshennik and D. K. Nadyozhin, Pis'ma Astron. 
Zh. 3, 353 (1977) [Sov. Astron. Lett. 3, 188 (1977)]. 

46. V. S. Imshennik and D. K. Nadyozhin, Usp. Fiz. Nauk 
156,576(1988). 

47. V. S. Imshennik and D. V. Popov, Pis'ma Astron. Zh. 
20, 620 (1994) [Astron. Lett. 20, 529(1994)]. 

48. V. S. Imshennik and D. K. Nadyozhin, Pis'ma Astron. 
Zh. 18, 195 (1992) [Sov. Astron. Lett. 18, 79(1992)]. 

49. V. S. Imshennik and D. K. Nadyozhin, Itogi Nauki 
Tekh., Ser. Astron. 21,63(1982). 

50. V. S. Imshennik and D. V. Popov, Preprint MPA 940 

(1996) . 

51. V. S. Imshennik and D. V. Popov, Pis'ma Astron. Zh. 
24, 251 (1998) [Astron. Lett. 24, 206(1998)]. 

52. V. S. Imshennik and D. V. Popov, Pis'ma Astron. Zh. 
28, 529 (2002) [Astron. Lett. 28, 465 (2002)]. 

53. L. N. Ivanova, V. S. Imshennik, and D. K. Nadyozhin, 
Nauchn. Inform. Astron. Soveta Akad. Nauk SSSR 
13, 3 (1969a). 



54. L. N. Ivanova et at., Proc. Int. Seminar on Physics 
of Neltrino and Neitrino Astrophysics (FIAN 
SSSR, Moscow, 1969), Vol. 2, p. 180 [in Russian]. 

55. E. Kolb et at., Phys. Rev. C 49, 1122(1994). 

56. L. D. Landau and E. M. Lifshits, Field Theory (Nau- 
ka, Moscow, 1973) [in Russian]. 

57. R. Mashuv, Prog. Part. Nucl. Phys. 40, 183 (1998). 

58. S. L. Mintz et al, Phys. Rev. C 40, 2458 (1989). 

59. D. K. Nadyozhin and N. V. Otroshchenko, Astron. Zh. 
263,366(1982). 

60. D. K. Nadyozhin, Astrophys. Space Sci. 49, 399 
(1977a). 

61. D. K. Nadyozhin, Astrophys. Space Sci. 51, 283 
(19776). 

62. D. K. Nadyozhin, Astrophys. Space Space Sci. 53, 
131 (1978). 

63. B. Paczynski, Annu. Rev. Astron. Astrophys. 9, 183 

(1971) . 

64. B. Paczynski and R. Sienkiewicz, Acta Astron. 22, 73 

(1972) . 

65. R. C. Peters, Phys. Rev. 136, B 1224 (1964). 

66. P. C. Peters and J. Mathews, Phys. Rev. 131, 435 
(1963). 

67. G. Pizzella, Nuovo Cimento B 102, 471 (1989). 

68. A. De Rujula, Proc. EPS Conf., Uppsala, 1987. 

69. A. De Rujula, Phys. Lett. B 193,514(1987). 

70. O. G. Ryazhskaya and V. G. RyasnyT, Pis'ma Zh. 
Eksp. Teor. Fiz. 56, 433 (1992) [JETP Lett. 56, 417 
(1992)]. 

71. O. G. Ryazhskaya and V. G. Ryasny, Nuovo Cimento 
A 106, 257(1993). 

72. J. L. Tassoul, Theory of Rotating Stars (Princeton 
Univ. Press, 1978). 

73. A. V. VoevodskiT, V. L. Dadykin, and O. G. Ryazh- 
skaya, Prib. Tekh. Eksp. 1, 85(1970). 

74. P. Yaranowski and A. Krolak, Astrophys. J. 394, 586 
(1992). 

75. J. R. Wilson et al., Ann. N.Y Acad. Sci. 470, 267 
(1986). 

76. Ya. B. Zel'dovich and O. Kh. Guseinov, Dokl. Akad. 
Nauk SSSR 162, 791 (1965). 

Translated by V. Astakhov 



